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ABSTRACT 

The  scope  of  this  thesis  is  to  establish  design  rules  for  the  shielded  form  of  Sus¬ 
pended  Substrate  Stripline  (SSL)  as  a  propagation  medium,  and  also  to  introduce  a 
practical  model  for  calculation  of  the  gap  dimensions  of  the  SSL. 

The  results  of  this  study  show  that  there  are  available  accurate  and  simple  design 
formulas  for  the  analysis  and  synthesis  calculations  of  SSL  transmission  line  parameters. 
It  is  also  shown  that  the  proposed  method  for  analysis  of  the  gap  discontinuity  can  be 
used  for  practical  applications  with  accurate  results. 
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I.  INTRODUCTION 


3  Planar  and  integrated  microwave  circuitry  for  application  in  radar  and  guidance 
systems  is  predominantly  based  on  the  microstrip  transmission  medium.  Due  to  high 
micrcfstrip  losses  with  increasing  frequency,  lower-loss  media  are  required  for  frequencies 
at  -band  (30GHz)  and  above.  These  higher  frequencies  are  coming  into  increased  use 
for  higher  target  resolution  and  improved  security. 

The  important  alternative  microwave  transmission  media  which  preserve  the  ad¬ 
vantages  of  planar  printed  circuit  forms  are  the  finline  and  shielded  suspended-substrate 
systems.  Of  those  two  media,  the  finline  has  had  extensive  application  and  theoretical 
investigation.  The  shielded  suspended-substrate  stripline  (SSL)  has  the  potential  for  very 
low  transmission  loss  since,  in  its  quasi-TEM  propagation  mode,  the  current  densities 
are  very  low  due  to  its  additional  air  gap,  relative  to  microstrip.  The  SSL  is  at  present  a 
developmental  system,  and  well-verified  design  rules  are  not  yet  available  for  it. 

A  number  of  empirical  formulas  for  the  computation  of  the  SSL  transmission  line 
parameters  have  been  published  in  the  literature,  but  the  validity  and  accuracy  of  these 
formulas  is  not  known.  At  present,  essentially  no  information  is  available  for  disconti¬ 
nuities  such  as  capacitive  gaps  or  angles  in  this  medium.  The  series  capacitive  gap  is  of 
particular  interest  in  that  a  useful  microwave  filter  can  be  built  using  resonator  sections 
coupled  by  scries  capacitive  gaps,  /tj — -  £\/vm/aM 

Of  the  candidate  theoretical  tools  Available  for  the  computation  of  electromagnetic 
field  problems:  the  method  of  moments  (Galerkin  method),  and  the  variational-solution 
method.  The  variational  method  constitutes  the  most  suitable  approach  for  obtaining 
the  wanted  theoretical  confirmation  of  SSL  TEM-mode  transmission  parameters.  Only 
a  limited  amount  of  preliminary  analytical  reduction  is  required  for  this  method  and 
well-developed  variational  formulas  are  available  in  the  literature  for  the  line  geometry 
and  the  discontinuities  of  interest. 

The  scope  of  this  thesis  is  to  compare  the  available  methods  for  analysis  of  SSL 
transmission  line  parameters  and  to  select  a  suitable  approach  for  practical  computa¬ 
tions,  and  also  to  examine  and  implement  a  practical  model  for  calculation  of  the  gap 
capacitance  in  the  SSL. 

The  verification  of  the  SSL  transmission  formulas  which  is  the  objecti.c  of  this  re¬ 
search  will  facilitate  the  application  of  this  low  loss  transmission  medium  in  radar  and 


millimeter  wave  system  design  and  contribute  to  the  available  knowledge  of  the  factors 
related  to  the  practical  utilization  of  the  system. 

The  thesis  begins  in  Chapter  2  with  the  examination  of  available  models  for  the 
analysis  of  transmission  line  parameters  of  shielded  SSL.  After  a  comparison  between 
the  models,  the  most  suitable  model  for  implementation  is  proposed. 

In  Chapter  3  the  calculation  of  static  capacitance  of  a  series  gap  in  shielded  SSL  is 
presented,  based  on  a  theoretical  method  described  in  (Ref.  4]. 

In  Chapter  4  the  final  SSL-filter  design  is  shown,  based  on  the  basic  filter  theory  and 
on  the  results  from  Chapters  2  and  3.  In  this  chapter  also,  the  designed  filter  is  evaluated 
by  using  the  Touchstone  CAD  Program  available  in  the  Microwave  Laboratory  of  Na¬ 
val  Postgraduate  School. 

Chapter  5  includes  the  conclusions  and  recommendations  of  this  study. 
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II.  MODELS  FOR  ANALYSIS  OF  SSL  TRANSMISSION  LINE 

PARAMETERS 


The  objective  of  this  chapter  is  to  examine  published  closed-form  design  formulas 
•  for  the  Suspended  Stripline  (SSL)  parameters,  and  to  select  accurate  and  practical  ex¬ 

pressions  for  implementation. 

Attention  in  this  work  is  restricted  to  TEM-mode  propagation  parameters.  This 
limitation  excludes  the  dispersive  (frequency-dependent)  effects  which  result  from 
hybrid-mode  wave  propagation,  which  occurs  in  mixed-dielectric  media.  This  simplifies 
the  analytical  task,  and  experience  has  shown  that  the  TEM  model  is  satisfactory  for 
microwave  frequencies  up  to  10  Ghz  and  is  frequently  extended  to  higher  frequencies. 

In  section  A  two  published  closed-form  design  models  for  shielded  SSL  are  exam¬ 
ined,  which  introduce  design  formulas  for  the  shielded  transmission  line  parameters.  The 
dominant  interest  in  this  investigation  is  the  shielded  form  of  suspended  striphne  typi¬ 
cally  operated  within  a  standard  waveguide  as  the  shielding  element.  The  advantages  of 
shielded  operation  for  microwave  transmissions  are  that  it  leads  to  reduced  radiation  and 
limited  interaction  with  neighboring  components.  In  section  B  two  available  fundamen- 
»  tal  methods  for  analysis  of  SSL  are  examined  which  are  based  on  electromagnetic  field 

theory.  These  are  used  to  test  the  formulas  cited  in  part  A.  In  section  C,  the  published 
/  empirical  design  formulas  for  unshielded  SSL  are  tested  and  the  results  are  compared 

with  the  formulas  for  shielded  SSL.  Finally,  in  section  D,  the  most  suitable  formulas  for 
the  foregoing  procedures  are  proposed  for  the  computation  of  SSL  transmission  line 
parameters. 


A.  CLOSED-FORM  DESIGN  EQUATIONS 
1.  Analysis  Calculations 

The  single  available  published  source  of  closed-form  expressions  for  shielded 
suspended  stripline  is  that  of  Ref.  1,  p.  693. 
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In  the  reference  closed-form  analysis  equations  for  shielded  SSL  are  developed, 
by  using  least-square  curve  fitting  to  numerical  results  of  the  finite-difference  methoJ. 


f - »  - » 


Figure  1.  Cross-Sectional  View  of  the  Transmission  Line  to  be  Analyzed. 


The  design  equations  of  Ref.  1  are  shown  below: 
The  characteristic  impedance  can  be  expressed  as 

Z 


(1) 


where: 

•  Z0  is  the  characteristic  impedance  of  SSL  of  identical  dimensions  and  completely 
filled  with  air,  and 

•  t,  is  the  effective  dielectric  constant  of  SSL. 

The  effective  dielectric  constant  is  given 


V^  =  [l+(£-£In-^-)ln 


where,  for  0  <  W  <  a/2, 

£  =  0.2077+  1.2177(4 ) -0.08364( -f ) 
b  b 


(2) 


(3) 


£=  0.0345 1  —  0. 103 1  ( j )  +0.0 1 742(  j ) 


(4) 
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or,  for  a/2  <  w  <  a, 


£  =  0.4640  +  0.9647(  4  )-0.2063(-f)  (5) 

b  b 

F=  -0.1424  +  0.301 7(  ^ )  -0.0241 1(  y )  (6) 

The  characteristic  impedance  is  given  : 


•  for  0  <  VV  <  a/2 


(7) 

where 

tj0  =  1207T 

(8) 

V -  -  1 .7866  -  0.2035(  -f )  +  0.4750(  •£ ) 

D  D 

(9) 

R  «  1 .0835  +  0. 1007(  y )  -  0.09457(  y ) 

(10) 

•  for  a/2  < 

W  <  a 

v^Zo  =  nQ[_V  +/?[  —  +  1.3930  +  0.6670  ln(  -y  +1.444)]-1] 

(ID 

w'here 

V=  —  0.6301  -  0.070S2(  4 )  +0.2470(  ■£ ) 

b  b 

(12) 

R  =  1.9492  +  0.1 553(~)  —0.51 23(-f ) 
b  b 

(13) 

where  W,  h,  a,  b  are  as  shown  in  Fig.  1. 

The  range  of  structural  and  substrate  constants,  over  which  the  developed 
equations  are  valid  is: 

1  <  a,b  <  2.5 
1  <  £,  <  4 
0.1  <  h/b  <  0.5 


2.  Synthesis  Calculations 

In  this  work  the  term  synthesis  is  used  to  denote  the  computation  of  trans¬ 
mission  system  dimensions  (Ref.  2:  p.  331)  to  achieve  specified  transmission  parameters, 
Z  and  i, ,  in  keeping  with  the  convention  ertablished  in  transmission  design.  The  com¬ 
putation  of  parameters  of  a  line  of  given  dimensions  is  called  the  analysis  process. 

In  the  method  of  Ref.  2,  a  set  of  simple  and  explicit  synthesis  equations  for  SSL 
was  developed  using  least-square  curve  fitting  to  numerical  results  obtained  using  Super 
Compact  for  SSL. 

For  specified  characteristic  impedance  Z  of  SSL,  the  strip  width  w  can  be  syn¬ 
thesized  using  the  following  formula 


’Y  -  [A  exp(  1.77245  -g- )  +B  exp(  -1.77245  -^§-  )]-'  (14) 


Figure  2.  Cross  Section  of  SSL. 


where  A,  B,  C,  D  have  different  values  corresponding  to  the  dielectric  constants  of  the 
substrate  as  follows: 

1.  For  c,  =  2.22 

A  =  0.0854(  -f  )"3  +0.290 1(  ~ )  +0.2342  -  ( )35(  -f  -0.4)  (15) 

b  b  Zb  b 

where 

S{  -7-  -0.4)  =  1  for  hi bZ  0.4  (16) 

b 

=  0  for  h/b  <  0.4 
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B  =  -  1 .2083 1  ln(  )  -0.3796  ln(  j )  +0. 1 370  ( 1 7) 

C=  1.446 
D  =  2.010 

2.  Fore,  =  2.80 

A  =  [0.09680  ln(y)  +0.61 30]/[0.3980  ln(  j )  +1.0276]  (18) 

B  =  1.038(  -f )  — 0.754(  -f )  -0.7380  ln(  ~ )  -0.4723  (19) 

b  b  b 

C=  1.440 
D  =  2.005 

3.  For  er  =  3.80 

A  =  [0.1 138  ln(  4- )  +  0.5775]/[0.4902  log(  -f )  +1.0238]  (20) 

0  D 

B=-  1.6943  ln(  \ )  - 1.544  -  ( f  -l)/(3.5  -|- )  (21) 

C=  1.685 
D  =  2.792 

The  synthesis  equations  of  SSL  are  valid  over  under  following  conditions: 


40£2  <>Z<>  1500 

(22) 

l<-f<2.5 

b 

(23) 

0.1  £-7-  £0.5 
b 

(24) 

3.  Implementation 

The  closed-form  equations  of  the  sections  above  were  evaluated  in  WR-28 
waveguide  with  z,f!-  2.22  and  substrate  height  h  =  0.5  mm. 

Figure.  3  shows  the  resulting  calculations  of  SSL  characteristic  impedance  Z0, 
for  stripwidths  W  from  1  to  7  mm. 

The  discontinuity  in  the  "Analysis  curve"  at  W = A/2  occurs  at  the  transition 
between  the  two  sets  of  design  equations  for  the  ranges  of  0  <  IV  <  —  and 

f <W<A- 

As  is  shown  in  the  same  figure  the  agreement  between  the  Analysis  and  Syn¬ 
thesis  calculations  is  within  1%  for  the  whole  range  of  W. 
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Figure  3.  Comparison  between  Analysis  and  Synthesis  Calculations. 
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B.  ELECTROMAGNETIC-FIELD  MODELS  FOR  SSL 

The  theoretical  models  used  in  this  comparison  were  a  variational  calculation  based 
on  Ref.  3,  p.  238  and  a  Quasi-static  solution  of  Laplace's  equation  based  on  Ref.  4  with 
boundary-point  matching  over  a  number  of  points  on  the  dielectric  surface  containing 
the  strip  conductor. 

The  variational  method  uses  Green's  functions  for  formulating  the  problem  and  a 
variational  principle  for  obtaining  line  capacitances. 

The  first  method  assumes  the  charge  density  across  the  strip  conductor  as  a  trial 
function  with  maxima  at  the  edges  of  the  strip,  while  the  second  method  specifies  only 
the  boundary  potentials,  a  consideration  which  simplifies  the  problem. 

1.  *  Variational  Calculation  Method 

In  this  method  a  various  trial  surface  charge  distributions  were  used  for  the 
charge  on  the  conducting  strip  [Ref.  3:  p.  238]. 


Figure  4.  Cross-Sectional  View  of  the  Transmission  Line  to  be  Analyzed. 
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The  static  capacitance  per  unit  length  of  a  TEM  transmission-system  config¬ 
uration  can  be  obtained  from  a  calculation  of  the  electrostatic  field  energy  WE  and  with 
use  of  the  expression: 

Q2 

we-Jc  (25) 

where  Q  is  the  charge  per  electrode  and  C  the  unknown  capacitance. 

The  energy  WE  can  be  computed  from 

KWJJp<Mv  (26) 

where  <1>  is  the  electrostatic  potential  and  p  is  the  charge  density  on  the  strip  conductor. 

With  the  z-coordinate  taken  along  the  axis  on  the  strip  transmission  line  the 
potential  O  in  the  transverse  plane  can  be  obtained  from  the  charge  distribution  by 
means  of  a  Green's  function 

x',y')dx'dy'  (27) 


Suitable  forms  of  Green's  function  for  the  SSL  geometry  are  available  in  the 
literature  [  Ref.  3:  p.  239],  and  for  reference  are  given  in  the  Appendix  B. 

The  charge  Q  in  Equation  (25)  is  calculated  from 

2  =  (28) 


where  p  is  the  charge  distribution  on  the  strip  conductor.  This  distribution  is  not  known 
a  priori,  but  for  the  strip  conductor  is  known  to  be  concentrated  at  both  edges,  relative 
to  that  at  the  center. 
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For  the  charge  density  the  following  form  of  the  trial  function  is  assumed: 


/M— £■[!+*  |-jjr(*-f  II1]  for  f— f-SJrsf  +  -f-  (29) 


=  0 


otherwise 


where  k  is  a  constant  to  determine  the  shape  of  f(x).  This  constant  can  be  chosen  such 
as  to  maximize  the  line  capacitance  C,  in  the  sense  of  the  variational  calculation. 

With  the  capacitance  C  per  unit  length  of  the  line  available,  the  transmission¬ 
line  parameters  then  may  be  calculated  as  follows 


Characteristic  impedance 


Z 


1 

3.V10VCQ 


(30) 


Transmission  line  wavelength  at  frequency  f: 


where  C0  is  the  value  of  capacitance  when  the  dielectric  layer  is  replaced  by  air.  The 
complete  calculation  from  [Ref.  3]  is  given  in  Appendix  A,  for  reference. 

Equations  (30)  and  (31)  were  calculated  for  WR-28  waveguide,  and  Fig.  6  shows 
the  resulting  calculations  of  the  SSL  impedance  Z0  for  c,  —  2.22,  substrate  height  h  =  0.5 
mm  and  linewidth  W  between  1  and  7  mm, 


11 


2.  Quasi-static  Solutions  of  Laplace's  Equation 

The  necessity  for  constructing  a  postulated  charge  distribution,  as  required  for 
the  variational  method,  can  be  avoided  by  treating  the  capacitance  calculation  as  a 
standard  boundary-value  problem  in  Laplace's  equation. 

With  this  approach  [Ref.  4]  ,  a  series  solution  is  written  for  each  of  the  three 
regions  in  the  line  cross  section  as  in  the  figure  below. 


Figure  5.  Cross  Section  of  SSL. 

If  a  potential  of  1  Volt  is  assigned  to  the  strip  conductor,  with  0  Volts  on  the 
outer  shield  boundary,  then  the  capacitance  can  be  obtained  as  the  integral  of  the 
normal  D-vector  over  a  bounding  surface  enclosing  the  center  conductor  where 

D-  —  cV<f>  (32) 

The  assumed  solution  in  region  3  is 

AT 

^3  -  Yj°n  sinh(W'  “  h3))  cos (V)  (33) 

n 

where  the  coefficients  D„  are  evaluated  at  N  selected  points  over  the  surfase  containing 
the  strip  conductor.  The  number  of  boundary  matching  points  N,  determines  the  accu¬ 
racy  of  the  solution  and  the  size  of  the  matrix  to  be  solved, 
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Solutions  of  the  type  of  equation  (33)  also  hold  in  regions  (1)  and  (2).  The  co¬ 
efficients  of  these  solutions  are  evaluated  by  enforcing  boundary'  conditions  at  the 
interfaces  between  regions.  The  design  equations  from  [Ref.  4]  are  given  in  Appendix  B 
for  reference.  These  equations  were  tested  in  WR-28  waveguide  and  Fig.  7  shows  the 
resulting  calculations  of  the  SSL  impedance  Z0  for  c,  =2.22,  h=0.5  mm  and  W  be¬ 
tween  1  and  7  mm. 

3.  Comparison  Between  Closed  Form  and  E/M  Field  Model  Calculations 

The  Analysis  and  Synthesis  closed-form  equations,  described  in  part  A,  were 
tested  against  the  more  fundamental  calculations  based  on  the  electromagnetic  field 
models,  described  in  part  B. 

Figure  8  shows  the  resulting  comparison  between  those  four  calculations  of  the 
SSL  impedance  Z„ ,  for  e,  =  2.22,  substrate  thickness  h=0.5  mm  and  stripwidth  W  be¬ 
tween  1  and  7  mm.  As  shown  in  Fig.  8,  the  agreement  is  within  1%  for  W  from  1  and 
3  mm,  and  8%  or  better  for  W  from  3  to  7,  with  the  greater  accuracy  at  the  more  fre¬ 
quently  used  line  widths. 

C.  COMPARISON  OF  DESIGN  EQUATIONS  FOR  SHIELDED  AND 
UNSHIELDED  SSL. 

Although  numerous  closed-form  expressions  exist  for  the  unshielded  form  of  SSL 
[Ref.  5:  p.  173],  [Ref.  6:  p.  453],  [Ref.  7:  p.  1429],  it  was  found  by  comparative  calculation 
of  the  two  cases  (shielded  and  unshielded)  that  the  presence  of  shielding  has  a  major  ef¬ 
fect  on  the  propagetion  characteristics  of  the  line  and  that  the  line  parameters  may  differ 
by  25%  to  40%  between  shielded  and  unshielded  lines  of  similar  dimensions.  Figure  9 
shows  the  resulting  comparison. 

Therefore  the  unshielded-line  design  equations  are  not  valid  for  the  case  of  Shielded 
SSL. 

D.  CONCLUSION. 

In  view  of  the  results  reported  above  we  can  conclude  that: 

•  The  analysis  formulae  [Ref.  1:  p.  693]  and  the  synthesis  formulae  [Ref.  2:  p.  331] 
introduce  accurate  and  simple  relations  for  calculation  of  SSL  transmission  line 
parameters. 

•  The  formulas  for  the  ushieldcd  SSL  are  non-valid  for  use  in  the  shielded  SSL  cal¬ 
culations. 
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Figure  7.  Characteristic  Impedance  Estimated  Using  the  Quasistatic  Method. 
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Figure  9.  Comparison  Between  Design  Equations  for  Shielded- Unshielded  SSL. 
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III.  DISCONTINUITY  GAP  CAPACITANCE 

A.  INTRODUCTION 

Since  the  development  of  microwave  integrated  circuits,  strip-type  transmission  lines 
have  been  widely  used  as  fundamental  structures. 

Although  a  great  deal  of  work  has  been  published  on  the  properties  of  the  suspended 
strip  transmission  lines,  the  theoretical  or  experimental  results  have  been  almost  entirely 
limited  to  the  impedance  and  the  phase  velocity. 

Since  published  results  for  the  discontinuity  structures  in  suspended  strip  trans¬ 
mission  lines,  such  as  the  gap  or  an  abruptly  ended  strip  conductor,  are  not  now  avail¬ 
able,  especially  in  the  case  of  shielded  SSL,  there  appears  to  be  a  need  for  quantitative 
characterization  of  the  discontinuity  structures,  either  in  term  of  equivalent-circuit 
models,  or  scattering-parameter  characterization. 

The  purpose  of  this  chapter  is  to  present  a  theoretical  method  for  calculation  of  the 
static  capacitance  of  a  series  gap  in  Shielded  Suspended  Stripline  based  on  Ref.  4  and  to 
show  the  procedures  used  for  its  implementation. 

B.  CALCULATION  OF  GAP  CAPACITANCE  IN  SHIELDED  SSL  IREF.  4J 

The  physical  gap  structure  in  the  strip  conductor  of  the  microstrip  transmission  line 
in  shown  in  Fig.  10. 


STRIP  UNE 


Figure  10.  Physical  Gap  Structure. 
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The  individual  gap  is  expected  to  have  a  series  capacitance  between  the  open  ends 
of  line,  plus  capacitances  to  ground  (shield)  due  to  the  fringing  Helds  at  the  break  in  the 
line.  These  capacitances  can  be  calculated  as  follows: 

By  placing  an  element  of  the  line  in  a  shield  box  with  1/2  gap  length  at  each  end, 
as  shown  in  Fig.  11,  the  gap  structure  can  be  represented  by  an  equivalent  n  -circuit  as 
shown  in  Fig.  12. 


Figure  12.  Equivalent  n  Circuit. 
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where 

C,  =  series  gap  capacitance  which  represents  the  electrostatic  capacitance  between 
the  two  open  ends  of  line  at  the  gap.  This  is  represented  here  by  capacitor  2 Ct  and  the 
image  in  the  end  wall  of  the  shield,  as  shown  in  Fig.  12. 

C0  =  line  capacitance  of  the  uniform  microstrip  transmission  line  with  length  1. 

C,  =  parasitic  capacitance  due  to  electrostatic  field  lines  extending  from  the  open 
end  to  ground. 

By  finding  a  way  to  calculate  the  total  capacitance  between  the  line  section  and  the 
enclosing  shield  box  for  the  two  cases: 

•  with  electric  walls  at  the  end  (giving  CE). 

•  with  magnetic  walls  at  the  end  (giving  CM). 

then  the  gap  capacitance  Ct  and  parasitic  capacitance  CP  can  be  calculated  from: 


(34) 


CP- 


(35) 


The  method  for  obtaining  electric  and  magnetic  walls  at  the  end  of  the  box  is  to 
assume  an  electric  E-ficld  within  the  shield  which  has  no  component  perpendicular  to  the 
end  walls  (magnetic-wall  case),  with  the  opposite  being  true  for  the  electric-wall  case. 
These  assumed  fields  are  written  in  the  form  of  finite-series  summations.  Each  term  in 
the  series  is  separately  a  solution  of  Laplace's  equation,  making  the  assumed  potentials 
valid  electrostatic  fields.  A  separate  series  is  assumed  for  each  of  the  three  regions  within 
the  shield:  upper  air  region,  the  dielectric  layer,  and  the  lower  region.  Each  of  these  series 
contains  constant  coefficients  with  each  term,  which  must  be  determined.  The  method 
of  determining  the  coefficients  is  to  require  matching  of  electrostatic  potentials  at  the 
upper  interface  (containing  the  conductor)  and  matching  of  normal  D-vectors  at  the 
lower  dielectric-air  interface. 

A  scries  is  assumed  solution  for  the  potential  in  each  of  the  three  regions  1,  2,  3  as 
shown  in  Fig.  13. 


<I>3  =  sinli (kq{v"  -  //3))  cos {kmx)  cos(A'„z)  (36) 

n  m 
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(37) 


<D2  -  sinh(A/)  +  Cmn  cosh(/c/)]  cos(Amx)  cos (knz) 

n  m 


<D,  =  Y^£jAmn  sinh(A^)  cos(*mjr)  cos {knz)  (38) 

n  m 


Figure  13.  Potential  Representation. 


where  Hn  is  the  air  height  of  region  n,  and  ki  ,  km,  and  kn  are  given  below. 

With  strip  line  at  a  potential  of  1  Volt  and  shield  at  0  Volts,  integrating  the  D  vector 
over  a  closed  surface  surrounding  the  strip,  the  capacitance  line-to-ground  is  given  by: 


C  = 


-eoJJV/Vdys 


(39) 


where  i=  1,2,3. 


The  surface  of  integration  is  broken  up  into  the  parts  covering  the  regions  1,2  and 
3  where  the  corresponding  d>,  are  used. 

The  coefficients  in  the  Equations  (36),  (37),  and  (38)  are  given  as: 


,  (2w-1)tt 

m  A 

(40) 

(2n—\)n  . 

- - 2j -  ( Electric  wall) 

(41) 

kn  —  p  ( Magnetic  wall) 

(42) 

II 

+ 

(43) 

where  D  is  the  length  of  the  shield  box  as  shown  in  Fig.  11. 
Using  the  following  boundary  conditions: 

•  Interface  (l)-(2) 


<ftl>, 

dy 


0>i  =0>2 


{aty  =  Hx)  -=  ir  —r  {aty'  =  0) 


•  Interface  (2)-(3) 


d>2  =  <1>3 


for  |  .v  |  < 


IL 

2 


0,=  1 


for 


d<X>2  d(P3 

~=  dy" 


(44) 

(45) 


(46) 

(47) 

(48) 


the  coefficients  Am ,  Dmn ,  C„„  can  be  evaluated  in  terms  of  Dmn  using  equations  (44)-(48) 
at  the  interfaces  between  layers.  Then  all  potentials  are  known  in  terms  of  the  unknown 
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To  solve  for  these  coefficients  we  evaluate  boundary  conditions  on  a  grid  of  selected 
points  (xft  z,)  over  the  interface  containing  the  strip  ,  as  shown  in  Fig.  14. 


w 

-D 

2 

l 

12 

i>  '  0/2 

Figure  14.  Grid  of  Points. 


As  more  points  are  included,  the  accuracy  of  the  solution  increases.  The  boundary 
conditions  lead  to  a  set  of  matrix  elements  to  solve  for  the  unknown  coefficient  Dv  after 
coefficients  Dm„  are  mapped  into  the  single-subscript  unknown  £>t ,  as  a  matrix  calcu¬ 
lation: 


DW-W  -■ TO 


where 

•  for  xf,  zf  on  the  metal  strip 

Bm  -  S3  cos(/cmx p)  cos {krfp) 
Wp—  1 

•  for  x,,zf  on  the  air-dielectric  boundary 


Bpq  -  kq 


cos {kmxp)  cos (k„yp) 


Wp- 0 


(49) 


(50) 

(51) 


(52) 

(53) 


where 


KN  -  cosh(A'?//A-) 
SlV=  sinh(y/lV) 

S13  =  sinh(*?(//,  +  //3)) 


(54) 


(55) 

(56) 


where  N=  1,2,3  and  the  HN  are  the  heights  of  the  respective  regions. 
The  capacitance  integral  produces  the  formula: 


CAP  —  AtQ 


K\SiK$  +  srA^5|3  4- 
A-,S2  -t-  A^-S 1 


(57) 


where 


D 


P=Dq  -TT-  sin(^m  -T- )  sin(fcm  -y ) 

Km*n  L  *• 


(58) 


C.  IMPLEMENTATION 

Using  the  previously  described  method,  a  Fortran  program  was  developed  for  cal¬ 
culation  of  the  gap  and  the  parasitic  capacitances  of  a  given  filter  structure.  This  pro¬ 
gram  is  given  in  Appendix  B. 

The  basic  steps  for  the  program  development  are: 

1.  Create  a  mesh  of  points  (*, ,  z;)  (as  in  Fig.  14)  on  the  (x,z)  plane  and  store  them  in 
array  [.v,z], 

2.  Compute  the  coefficients  k„  (equation  41)  for  n=  odd  and  store  them  in  an  array 
[£„(«,)],  in  order  to  use  them  for  the  electric  capacitance  calculations.  Compute 
kn  for  n  =  even  (equation  42)  and  store  them  in  array  [£„(«:)]  to  use  for  the  mag¬ 
netic  capacitance  calculations, 

3.  The  array  of  points  xi}  z ,  has  dimensions  7mix  by  Jm„  .  It  is  necessary  to  convert  this 
to  one-dimensional,  string  by  using: 

P  ~  AnaxO  “  1)  +  i  for  the  ( ij )tb  member  (59) 

Similarly  the  double  subscripted  Dnn  has  to  be  converted  to  a  single-subscript 
quantity  £>, ,  to  enable  solution  of  the  matrix.  This  is  done  by  using  : 

<7  =  '»max("  -!)  +  »»  (6°) 
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4.  Compute  the  coefficients  km  (equation  40)  for  m=  1  to  mm„  and  store  them  in 

array  [AJ. 

5.  Compute  the  coefficients  A,  (equation  43)  for  A„(W|)  and  calculate  the  arrays 
[Z?(p,<7)]  and  [  IF(p)]  ,  on  the  metal  strip  and  on  the  air-dielectric  boundary. 

6.  Use  Gaussian  elimination  to  solve  for  the  unknown  coefficients  Z)f  as  : 

-[*„]-■  [«g. 

7.  Calculate  the  electric-wall  capacitance  using  the  equation  (57). 

8.  Compute  the  A,  (equation  43)  for  An(«2),  calculate  the  arrays  [£(/>,<?)]  and 
[fF(p)]  ,  and  solve  to  compute  the  Dr 

9.  Calculate  the  magnetic  capacitance  using  the  equation  (57). 

10.  Calculate  the  gap  capacitance  using  equation  (34). 

11.  Calculate  the  parasitic  capacitance  using  equation  (35). 

The  Fortran  program  given  in  Appendix  B  was  implemented  on  the  Naval  Post¬ 
graduate  School's  IBM  3033  computer,  for  maximum  number  of  points  (allowed  by  the 
available  storage  capacity),  7mlx  =  23,  Jmtx  =  23  and  also  =  23  and  -  23  . 

Because  of  the  limitation  of  the  number  of  points  available,  it  was  necessary  to  ad¬ 
just  the  distribution  of  the  points  in  the  (x,z)  plane,  in  order  to  have  more  points  in  the 
area  of  the  gap-strip  boundary.  This  was  done  by  using  the  formulas: 

for  (61) 

=  +  for  J>Jf  (62) 

where  Jf  is  the  number  of  points  choosen  on  the  strip.  For  example,  for  Jmn  =  23  and 
for  =  10  ,  there  are  10  points  on  the  strip  and  13  points  outside  the  strip. 

Also  it  was  found  that,  if  the  distribution  of  the  points  is  not  correct,  the  Fortran 
program  gives  negative  values  of  the  electric  (C£)  or  magnetic  ( CM )  capacitances. 

In  this  event,  the  following  steps  are  taken: 

•  In  case  of  negative  electric  capacitance,  the  number  of  points  on  the  strip  (  JF )  must 
be  increased. 

•  In  case  of  negative  magnetic  capacitance,  the  number  of  points  outside  the  strip 
must  be  increased  (  J,  to  be  decreased). 

Implementing  the  Fortran  program  with  shield  dimension  D=  12  mm.  number  of 
points  of  23x23,  and  using  the  dimensions  of  the  waveguide  of  interest  (WR-28),  the  gap 
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capacitance  is  calculated  for  gap  dimensions  between  0.2  mm  to  0.8  mm,  as  shown  in 
Fig.  15. 


Figure  15.  Gap  Capacitance  Calculations. 
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D.  COMPARISON  OF  THE  RESULTS 


In  order  to  check  the  accuracy  of  the  results  given  from  Fortran  program  (Appen¬ 
dix  B),  the  capacitances  of  the  SSL  structure  shown  in  Fig.  16  were  determined. 


The  gap  capacitances  of  this  structure  are  given  in  [Ref.  10:  p.75]  calculated  by  using 
a  Fortran  program  based  on  the  work  of  Smith  (Ref.  11:  p.  424 j  and  are  given  by: 

For  Gap* 0.107  mm  Capacitance* 0.0210  pF 
For  Gap  *  0.307  mm  Capacitance  *  0.00634  pF 

Using  the  Fortran  program  given  in  Appendix  B,  with  dimension  of  the  shield  box 
D*  10  mm  (Fig.  11)  the  calculated  values  are: 

For  Gap* 0.107  mm  Capacitance* 0.0320  pF 
For  Gap* 0.307  mm  Capacitance* 0.00580  pF 

Thus  as  can  be  seen  the  calculated  gap  capacitances  are  in  reasonable  agreement 
with  the  values  of  those  in  [Ref.  10:  p.75J,  which  were  also  approximate  calculations 
from  a  static  model. 

E.  CONCLUSION 

Considering  all  of  the  above,  we  can  conclude  that  the  model  given  in  (Ref.  4J  for 
the  calculation  of  the  gap  capacitance  in  shielded  SSL  structures  is  relatively  easy  to 
implement  and  gives  accurate  results. 
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IV.  SHIELDED  SUSPENDED  STRIPLINE  FILTER  DESIGN 

A.  INTRODUCTION  < 

Suspended  and  inverted  microstrip  lines  are  among  the  principal  transmission  media 
used  in  the  upper  microwave  and  lower  mm- wave  bands.  A  characteristic  aspect  of  these  i 

lines  is  that  the  presence  of  an  air  gap  between  the  substrate  and  the  ground  plane  typ¬ 
ically  reduces  the  effects  of  dispersion  on  the  propagation  constant  to  an  extent  that  the 
quasi-static  results  remain  useful  even  at  very  high  frequencies. 

The  suspended  strip  line  (SSL)  is  a  modified  version  of  the  microstrip  line.  Com¬ 
pared  with  the  normal  microstrip  line,  it  has  some  attractive  features,  such  as  lower  at¬ 
tenuation  and  larger  tolerance  of  fabrication  variances. 

B.  MICROWAVE  FILTER  THEORY 

The  network  synthesis  methods  for  filter  design,  generally  start  out  by  specifying  a 
transfer  function,  as  a  function  of  complex  frequency  p. 

The  TchebvschefT  and  maximally  fiat  transfer  functions  are  often  used  for  filter  ap¬ 
plications.  In  this  work  only  the  TchebvschefT  design  is  considered,  since  is  preferred 
because  of  its  sharp  cutoff.  The  attenuation  characteristic  of  the  Tchebyscheff,  or  "equal 
ripple"  is  shown  in  Fig.  1 7.  The  attenuation  <5m  is  the  maximum  db  attenuation  in  the 
pass-band,  while  cu,  is  the  equal-ripple  band  edge. 

Attenuation  characteristics  of  this  may  be  specified  mathematically  as  [Ref.  9:  Ap¬ 
pendix  C]. 


<5  =  10  log(l  +  c  cos2 {ncos~\ -{%-)))  <y£<Uj 

U/] 

nER  ~  10  log(l  +  t  cosh2(«  cosh“'(  — ■ )))  <u  S:  w, 

2 

^ AR  w  ,  (Pm~l) 

c  =  (antilog(  — ))— 1=  4p  - 


(63) 

(64) 

(65) 
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where:  8  “insertion  loss,  nER  “isolation,  and  n*  number  of  reactive  elements  in  the 
circuit.  The  basic  prototype  circuit  is  shown  in  Fig.  17. 


Figure  17.  Prototype  Circuit  (Ref.2,  Appendix  CJ. 


The  element  values  g0»  &»••«  Sm  &.+t  for  a  low-pass  prototype  must  alternate  between 
shunt  capacitor  and  series  inductor.  The  structure  can  begin  with  either  inductor  or 
capacitor.  The  index  associated  with  each  element  is  given  by  k.  Then  the  following 
equations  define  the  element  values  [Ref.  1:  p.  97j,  [Ref.  2:  Appendix  C). 


WjQ 

Y  =  %k 

(66) 

<y,4 

z0  ~gk 

(67) 

„  2(11 

8 1  Y 

(68) 

Aak-\ak 
gk  b/c-lSk-l 

(69) 

«k  ~  «n(  2n  n) 

(70) 
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(71) 


bk-  y2  +  sin2(-“*) 

V  -  sinh(  ~  ]n{coth  -p-jy ))  (72) 

<5M-20/og(-^li)-(p-l)2  (73) 


where  pm  is  the  max  VSWR. 

C.  TRANSFORMATIONS 
1.  High-Pass  Prototype 


The  low-pass  prototype  is  converted  to  a  high-pass  prototype  by  using  series 
capacitors  and  shunt  inductors  [Ref.  9:  Appendix  C]  as  shown  in  Fig.  18. 


Figure  18.  High-Pass  Prototype  [Ref.  9:  Appendix  C). 


The  element  values  are  to  be  found  using  the  following  equations 

f^iQ  _  1 

Vo  8k 

_  1 

Z0  8k 


(74) 

(75) 
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Equations  (63)  and  (67)  are  the  same  except  that  attenuation  <5  is  above  fre¬ 
quency  w,  and  isolation  n  is  used  below  frequency  to,. 

2.  Band-Pass  Prototype 

The  low-pass  prototype  filter  is  converted  to  a  band-pass  filter  [Ref.  9:  Appendix 
C]  simply  by  resonating  all  of  the  elements  at  the  center  (geometric  mean)  frequency  of 
the  pass-band. 

First  the  ripple,  bandwidth  and  number  of  elements  are  selected,  to  give  the 
maximum  insertion  loss  and  bandwidth  desired  in  the  pass-band.  This  produces  a  low- 
pass  prototype  as  in  Fig.  17. 

Then  given  upper  frequency  a)„  and  lower  frequency  <oL  the  center  frequency 
w0  is  calculated 

wo  “  (76) 

Then  inductors  are  selected  to  parallel  resonate  with  the  shunt  capacitors  as 

L, - J-  C2 - e,c.  (77) 

The  resonating  elements  simply  shift  the  response  up  to  a  higher  frequency. 
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D.  CAPACITIVE-GAP-COUPLED  TRANSMISSION  LINE  FILTERS 

The  SSL  filter  considered  here  is  a  capacitive-gap-coupled  transmission  line  filter 
(Ref.  8:  p.  441]. 

Figure  19  shows  the  structure  of  a  coupled-resonator  filter  consisting  of 
transmission-line  resonators  which  are  approximately  a  half  wavelength  long  at  the 
midband  frequency  g)0  and  which  have  series  capacitance  coupling  between  resonators. 


Figure  19.  Capacitive-Gap  Coupled  Filters  (Ref.  1:  p.  441]. 


In  this  case  the  inverters  are  in  the  form  of  series  capacitive  types.  These  inveiters 
tend  to  reflect  high  impedance  levels  to  the  ends  of  the  half-wavelength  resonators  and 
it  can  be  shown  that  this  causes  the  resonators  to  exhibit  a  shunt-type  resonance  (Ref. 
8:  p.  440].  Thus  the  filters  under  consideration  operate  like  the  shunt-resonator  type  of 
filter  whose  general  design  equations  (Ref.  8:  p.  440]  are: 


4>.i  _  In  o)  ~~ 
>o  V  2  *oSi®'i 


(78) 


JJJ+\  nIV  1 

>0  2w'l  Jspgj+i 


for  j  =  1  to  n  —  1 
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where 

g0,  :  Tchebyscheflf  coefficients  [Ref.  8:  p.  100]. 

W  :  the  fractional  bandwidth 

i  :  admittance  inverter  parameters 
K„  :  characteristic  admittance  of  the  line 

Assuming  the  capacitive  gaps  function  as  ideal  series-capacitance  discontinuities  of 
susceptance  B}JJ.X  we  have: 


BJj+ 1 

>o 


JJJ+ 1 
>o 


l-(- 


JJ+1 
V  ' 


(81) 


1  _i  2 Bi—]j  , 

6,  =  n  -  -T-  [  tan  ( — rr-^-+tan“( — ~ — )]  in  radians  (82) 

i  *o  h 

where  BtJ.x  and  0,  are  evaluated  at 

Using  the  above  susceptance  values  the  gap-capacitances  are  derived  using  the  fol¬ 
lowing  equation: 


<83> 

E.  SUSPENDED  STRIPLINE  FILTER  DESIGN 

The  filter  design  parameters  selected  for  this  design  are 

•  Center  frequency  =  33.25  Ghz 

•  Bandwidth  =  2.6  Ghz 

•  Percentage  Bandwidth  =  8% 

•  Number  of  poles  =  3  and  5 

•  Ripple  =  0.2  db 

The  design  of  the  filter  consists  of  two  major  parts: 

•  Calculation  of  the  transmission  line  parameters  of  the  filter  line  elements. 

•  Calculation  of  the  gaps  and  the  resonance  lengths  of  the  filter. 
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1.  Calculation  of  the  Transmission  Line  Parameters  of  the  Shielded  SSL 

For  the  waveguide  of  interest  WR-28  (A  =7. 11  nun,  B  =  3.56  nun),  with 
e,  =  2.22  ,  h  =  0.5  nun  and  Z=  60  £2  ,  using  the  equations  (2)  and  (14)  we  calculate: 

W  =  3.33  nun 

c,  =  1.093769  where  e,  =  effective  dielectric  constant.  So  finally  the  filter  struc¬ 
ture  is  as  shown  in  cross  section  in  Fig.  20. 


Figure  20.  SSL  Filter  Cross-Sectional  View. 


2.  Calculation  of  the  Gaps  and  Resonance  Lengths 
a.  Three  Element  Filler 

For  N=  3  and  r=0.2db  and  using  the  coelficients  from  [Ref.  8:  p.  100] 
g,=g3=  1.2275,  g:  =  1.1525,  g4=  1 

•  Using  the  design  equations  (78)  through  (83) ,  with  the  following  specifica- 

i 

tions: 

•  Bandwidth  BW  =  0.08 

•  Characteristic  impedance  Z  =  60  Q 

•  Center  Frequency  =  33.25  Ghz 

•  Effective  dielectric  constant  k=  1.093769 

we  calculate: 
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•  Gap  capacitances 

C0!  =  0.28436  pF 
C„  =  .08523  pF 
C2  3  =  .08523  pF 
CM  =  0.28436  pF 

•  The  length  of  each  element  of  the  filter  in  radians 

0,  =  2.7266600  rad 
02  =  2.9310686  rad 
6 j  =  2.7266600  r<w/ 

•  The  actual  length  of  each  element  in  mm 

L,  =  3.772  mm 
Lj  =  4.055  mm 
Lj  —  3.772  mm 

•  Using  Fig.  15  in  Chapter  3  we  calculate  the  required  gaps  for  the  previously  speci¬ 
fied  gap  capacitances. 

A01  =  0.70  mm 
Au  =  0.75  mm 
Ay  =  0.70  mm 

•  Using  the  Fortran  program  in  Appendix  B  we  calculate  the  parasitic  capacitances 
as  follows: 

C„,a\  —  0.44  10-“  F 
C,or\i  —  0.20  10"16  F 
C =  0.44  10'1S  F 

b.  Five  Element  Filter 

For  N  =  5  and  r=0.2  db  we  pick  the  coefficients  from  [Ref.  8:  p.  100] 
ft 1.3394,  ft -ft -  1.3370,  ft  =  2.1660,  ft  =  1.3394,  ft  =1.0 


Using  the  design  equations  (78)  through  (83)  we  calculate: 

•  Gap  capacitances 

C0)  =  0.2696  pF 
C,2  =  0.0755  pF 
CJ3  =  0.0592  pF 
C34  =  0.0592  pF 
C«  =  0.0755  pF 
C56  =  0.2696  pF 

•  The  lengths  of  each  element  in  rad 

0,  =  2.7507  rad 
02  =  2.9742  rad 
<?3  =  2.9941  rad 
04  =  2.9742  rad 
6S  =  2.7507  rad 

•  The  actual  length  of  each  element  in  mm 

L,  =  3.776  nun 
L:  =  4.083  mm 
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^3  =  4.1 11  mm 
Lt  =  4.083  mm 
L$  =  3.776  mm 

•  Using  Fig.  15  in  Chapter  3  we  calculate  the  required  gaps  for  the  previously  spec¬ 
ified  gap  capacitances: 

A01  =  0.70  mm 
A„  =  0.75  mm 
AJ3  =  0.80  mm 
A34  =  0.80  mm 
Au  =  0.75  mm 
A<6  =  0.70  mm 

•  Using  the  Fortran  program  given  in  Appendix  B  we  calculate  the  parasitic 
capacitances  for  each  case  as  follows: 

C^-  0.440  10-“  F 
C,„12  =  0.830  10-“  F 
C,„i3  =  0.166  10-“  F 
Cpcn  =  0.166  10-“  F 
C„M  =  0.830  10-“  F 
C,„56  =  0.440  10-“  F 

Considering  all  the  above  calculations  the  complete  filter  sructures  are  as 
shown  in  Fig.  21  and  Fig.  22. 
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Figure  21.  Complete  3  Element  Filter  Structure. 


o.7mm  0.75mm  0.8mm  0.8mm  0.75mm  0.7mm 

“*l  1*1*1  U-"*!  1*1*1  U_"*l  1*1*1  U_ 


Figure  22.  Complete  5  Element  Filter  Structure. 
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F.  FILTER  EVALUATION 

Using  the  previously  calculated  values,  the  two  Filters  were  characterized  using  the 
Touchtone  EEsof's  computer-aided  microwave  simulation  and  optimization  program 
(available  from  the  EEsof  Corporation,  Westlake  Village  California).  The  models  used 
in  the  Touchstone  programs  are  shown  in  the  next  Figs.  23  and  24. 


Figure  23.  3  Element  Filter  Structure. 


Figure  24.  S  Element  Filter  Structure. 
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From  the  TOUCHSTONE  calculation  the  frequency  responses  of  the  two  filters  are 
plotted  as  shown  in  Figs.  25  and  26.  As  shown  in  those  two  figures  the  filter  bandwidth 
and  the  center  frequencies  are  in  agreement  with  the  design  values. 

The  effect  of  the  parasitic  capacitances,  because  of  their  very  small  values,  was 
negligible  in  the  frequency  response  of  the  filter. 
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EEsof  -  Touchstone  -  Thu  Jun  01  13: 29: 04  1989  ■ 


□  03[  Sll]  +  D3l  S21] 
FILTER  FILTER 


Figure  25.  Frequency  Response  of  3  Element  Filter. 
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Figure  26.  Frequency  Response  of  5  Element  Filter. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  purpose  of  this  research  was  to  establish  design  rules  for  the  shielded  form  of 
Suspended  Substrate  Stripline  as  a  propagation  medium,  and  also  to  establish  a  practical 
model  for  calculating  the  gap  capacitance  of  the  filter  structure. 

The  objective  of  this  research  was  accomplished.  The  final  results  are  the  following: 

1.  For  the  calculation  of  the  Shielded  SSL  transmission  line  parameters,  the  Analysis 
[Ref.  1:  p.693]  and  Synthesis  [Ref.  2:  p.331]  design  equations  can  be  used  because 
these  relations  are  useful  and  accurate. 

2.  The  method  given  in  [Ref.  4]  for  the  calculation  of  the  static  capacitance  of  a  series 
gap  in  shielded  SSL  is  relatively  easy  to  implement  and  gives  accurate  results. 

B.  RECOMMENDATIONS 

Based  on  the  results  of  this  study  the  following  recommendations  are  made: 

1.  The  dispersive  (frequency-dependent)  effects  must  be  addressed,  especially  when 
the  Suspended  Stripline  filter  is  to  be  used  for  microwave  frequencies  higher  than 
10  Ghz.  For  this  purpose  a  spectral-domain  analysis  of  the  propagation  in  this  line 
system  should  be  made,  based  on  wave  potentials  which  are  solutions  of  the 
Helmholtz  equation  in  the  guiding  region. 

2.  The  gap-capacitance  calculations  should  similarly  be  carried  out  in  a  frequency- 
dependent  spectral  domain  form.  Ultimately,  other  discontinuities,  including  the 
line  open  end  and  corner  (bend)  should  be  investigated. 
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APPENDIX  A.  DESIGN  EQUATIONS 

A.  VARIATIONAL  CALCULATION  METHOD 

This  is  a  general  method  (Ref.  3:  p.  238]  for  analyzing  the  transmission  line  charac¬ 
teristics  of  SSL,  with  rectangular  outer  conducting  boundaries.  This  method  uses  a 
Green's  function  for  formulating  the  problem  in  a  variational  approach  for  obtaining 
practical  solutions. 


Figure  27.  Cross  Section  of  SSL. 

The  characteristic  impedance  Z„  of  the  line  is  given  by 


cco 


(85) 


The  waveguide  wavelength  /  is 


and  the  wave  propagation  velocity  V  is  given: 

v  -  I 'oVlF  (*7) 


where: 

F0  “velocity  of  light  in  free  vacuum  space 
=  wavelength  in  free  vacuum  space 
C0  =  line  capacitance  for  the  case  c,  =  c2  =  £3  =  £0 
The  line  capacitance  is  given  : 


C— tt - - -  (88) 

j  J/WGCx,  *!  +  *2  |  +  h2]f[x0)dxdx0 

where  fl;x)  is  the  charge  density  distribution  across  the  strip  conductor  and  G(x,y)  is  the 
Green's  function  which  is  regarded  as  the  potential  at  (x,y),  due  to  a  unit  charge  in  an 
infinitely  small  volume  at  (^r0»  J’o)- 

The  charge  density  of  thin  strip  conductor  is  known  to  be  large  at  both  edges,  rela¬ 
tive  to  that  at  the  center.  For  this  reason  the  following  form  of  the  trial  function  is  as¬ 
sumed  in  this  method: 

/W— +  )|3]  for  f — y-SJTSf  +  -J-  (89) 


=  0  otherwise 

The  constant  K  is  chosen  such  that  to  maximize  the  value  of  the  capacitance  C. 
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The  line  capacitance  is  evaluated  by  applying  the  equation  (89)  in  equation  (88)  and 
is  given  by: 


where: 


OO 


K—w(Yu+kV. a,) 


e-1+f 


nnlV  \ 

v  _  2 IV  .  nn  >.  sm(  2 a  }  .  3 
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*  w^/j]  nr.lu  •  .  ,  ,  nnh,  nnh7 

'In  =  ci  cosh(  -j- )  smh(  — )  +  c2  smn(  — ~ )  cosh(  ) 


»  nnh3  .  nnh3  , 

=  e2c„  sinh(  — )  +  £3n„  cosh(  — — ) 

y  »  ,  .  n7th}  nnh2  ,  .  .  ,  ,  «7t/i,  .  nnlh 

£n  =  £]  cosh(  — -  ,  cosh(  — — )  +  £2  sinh(  ~5- )  sinh(  —^L ) 
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The  value  of  the  constant  K  is  given  by: 


Yjmn-Vu)V^, 

K - Jg -  (99) 

nm  1 

Implementation:  The  described  method  was  applied  to  a  specific  waveguide  of  interest 
(WR-28)  with  dimensions  a=7.11  mm  and  b=3.56  mm,  substrate  h=0.5  mm  and  t, 
■  2.22.  Equation  (99)  from  Ref.  3,  p.  238  was  not  used  because  it  was  found  that  it  gave 
unreasonable  values  (some  negative)  of  K.  Instead  of  equation  (99)  the  K  was  chosen 
such  as  to  maximize  the  line  capacitance  C  (in  this  case  K  =  3). 

B.  ALTERNATE  ASSUMPTION  FOR  CHARGE  DISTRIBUTION 

This  version  [Ref.  4]  of  the  variational  calculation  method  assumes  the  charge  den¬ 
sity  function  across  the  strip  to  be  simply  a  step  function  as  shown  on  Fig.  28. 


.W 


Figure  28.  Charge  Density  Distribution. 


By  proper  choice  of  the  pulse-function  height  and  width  at  the  edges,  a  good  repre¬ 
sentation  of  the  typical  charge  singularity  which  occurs  at  the  edges  of  strip  conductors 
was  obtained.  The  ratio  of  the  width  of  charge  peaks  to  linewidth  is  defined  as: 


o 


(100) 
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The  line  capacitance  is  given  by: 


(107) 


JW*  =  W\_l+2a(Q~l)2  (108) 

Design  approach:  Considering  all  the  above  the  design  approach  is 


q  num 

denoni 

(109) 

OO 

Denom  -  ^£„(y)(/«;)2 

M=1 

(HO) 

(-^)l(2-l)cos(^)-cos(^f )  + 

+2(  cos( )  _  C0S(  }), 

(111) 

Mim- (Mil  +2<,(g- l)]j! 

(112) 

E„(y)  =  C„sinh(-f-  (*->•)) 

(113) 

D  _  2,,“ 

"  ^(At 

(114) 

»/„  as  given  in  Equation  (96) 
cn  as  given  in  Equation  (98) 

A„  as  given  in  Equation  (97) 

C  =  — 

C  2 

(115) 

II 

O 

l 

(116) 
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E=D  +  oW 

(117) 

G=C  +  - f- 

(118) 

F=  G  —  alV 

(119) 

The  values  of  Q  and  a  must  be  chosen  in  such  a  way  to  maximize  the  calculated  , 
capacitance  value. 

Implementation:  This  method  was  applied  to  a  specific  waveguide  of  interest 
(WR-2S)  with  dimensions  a  =  7. 1 1  mm,  b  =  3.56  mm,  h2  =  0.5  mm,  and  cr  -  2.22. 

The  parameters  Q  and  a,  were  chosen  in  such  a  way  to  maximize  the  calculated 
capacitance  and  for  this  specific  case  had  the  values  Q=  22  and  a  =0.01. 
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APPENDIX  B 


FORTRAN  PROGRAM 


C  THIS  PROGRAM  CALCULATES  THE  GAP  AND  THE  PARASITIC  CAPACITANCES 
C  OF  A  GIVEN  FILTER  STRUCTURE,  USING  THE  METHOD  DESCRIBED  IN 
C  CHAPTER  3. 

c 

DIMENTION  X(40) ,Y(40) ,KM(40) ,KN(60) ,B(485,485) ,D(485) ,W1(485) 
real  A,B1,H1,H2,H3,D1,G,ER,W,PI,F1,F2,K1,K2,K3,S1,S2>S3,S13 
REAL  G1,G2,G3,HH1,T,KQ,S6,C1M,C2M,C3M, C4M , NCAPM , CAPM 
REAL  C1E , C2E , C3E , C4E , NCAPE , CAPE , CO 

INTEGER  IMAX,MMAX,NMAX, JMAX, I, J,N1,N2,M2,M,N,I2, J2,I3,J3,L 
INTEGER  P,Q,K,V,N3,M3,Q1,Q2,N4,N5,N6,M5,N8,I4, J4,I5, J5,IP, JP 
c  A=width  of  waveguide 

c  Bl=length  of  waveguide 

c  H2=substrate  height 

c  W=stripwidth 

c  ER=dielectric  constant 

c  CO=capacitance  with  the  dielectric  layer  replaced  by  air 
c  Dl=length  of  the  shield  box 

c  G=gap  dimentions 

c  Imax=number  of  points  in  z-axis 

c  Jmax=number  of  points  in  x-axis 

c  Mmax,Nmax=summation  terms 

C  ******* INPUTS  ********* 

A=7.  11/iOOO 
Bl=3. 56/1000 
H2=0. 5/1000 
W=3.  33/1000 
ER=2.  22 
C 

C0= 

Dl=20/1000 
G=l.  0/1000 
IMAX=16 

C  ***  BALANCE  BETWEEN  THE  POINTS  *** 

C  IP=10 
JP=8 

C  ********************************** 

JMAX=IMAX 
MMAX=JMAX 
NMAX=IMAX 
PI=3.  1415926 
Hl=(Bl-H2)/2 
H3=H1 
Fl=W/2 
F2=(Dl-G)/2 
C 

V=IMAX*JMAX 

C  ***  Balance  in  x-axis  *** 

DO  50  1=1 , IMAX 
C  IF  (I  .LE.  IP)  THEN 
C  X( I )=( ( A-W) /( 2*IP) )*I 

C  ELSE  IF  (I  .GT,  IP)  THEN 
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C  X(I)=W/2+((A-W)/(2*IMAX))*I 

C  ENDIF 

X(I)=A*I/C2*IMAX) 
c  ***  Balance  in  z-axis  *** 

DO  45  J=1,JMAX 

IF  (J  .LE.  JP)  THEN  t 

Y(J)=((D1-G)/(2*JP))*J 
ELSEIF  (J  .GT.  JP)  THEN 
Y(J)=((Dl-G)/2)+(G/C2*JMAX))*J 
ENDIF 

45  CONTINUE 
50  CONTINUE 

C  **Calculation  of  coefficiets  Km  and  Kn  for  n=even  and  odd** 

DO  100  M=1,MMAX 
KM(M)=(2*M-1)*PI/A 
DO  80  N=1,NMAX 
N1=2*N-1 
N2=2*N 

KN(N1)=(2*N1-1)*PI/D1 
KN(N2)=2*N2*PI/D1 
80  CONTINUE 
100  CONTINUE 

C  **Calculation  of  coefficients  B(p,q)  and  W(p)** 

DO  300  12=1, IMAX 
DO  250  J2=l, JMAX 
DO  200  M2=l , MMAX 
DO  150  N2=1,NMAX 
P= JMAX* (12-1) +J 2 
Q=NMAX* ( M2 - 1 ) +N2 
KQ=( KM( M2 ) **2+KN( N2 )**2 ) **0 .  5 

IF  (X( 12)  .LE.  FI  .AND.  Y(J2)  .  LE.  F2)  THEN  . 

S3=SINH(KQ*H3) 

B(P,Q)=S3*COS(KM(M2)*X(I2))*COS(KN(N2)*Y(J2)) 

W1(P)=-1.0  * 

ELSE 

K1=C0SH(KQ*H1) 

K2=COSH( KQ*H2 ) 

K3=COSH(KQ*H3) 

S1=SINH(KQ*H1) 

S2=SINH(KQ*H2) 

S3=SINH( KQ*H3 ) 

S13=SINH(KQ*(H1+H3)) 

G1=ER**2*S 1*S2*S3 

G2=K1*S2*K3+ER*K2*S13+G1 

G3=COS( KM( M2 ) *XC 1 2 ) )*C0S  C  KN( N2 ) *Y( J2 ) ) 

G4=K 1*S  2+ER*K 2*S 1 
BCP,Q)=KQ*G2*G3/G4 
W1CP)=0. 0 
ENDIF 

C  WRITEC2, 130)  P,Q,B(P,Q) ,W1(P) 

C  130  FORMAT(I3,1X,I3,1X,F20.  15,1X,F20. 15) 

150  CONTINUE 

200  CONTINUE  4 

250  CONTINUE 
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300  CONTINUE 

C  **Solve  the  equation  (49)  in  order  to  compute  the  Dq** 
DO  655  13=1, V-l 
B1=ABS(B( 13,13)) 

B1=ABS(B1) 

L=I3 

14=13+1 

DO  350  J3  =  14, V 
IF  (ABS(B( J3,I3))  .LT.  Bl)  THEN 
GO  TO  350 

ELSEIF  (ABS(B( J3,I3))  .GE.  Bl)  THEN 
B1=ABS(B( J3,I3)) 

L=J3 
END  IF 
C 

350  CONTINUE 

IF  (Bl  .EQ.  0.0)  THEN 
GO  TO  850 

ELSEIF  (L  .EQ.  13)  THEN 
GO  TO  500 

ELSEIF  (L  .NE.  13)  THEN 
GO  TO  400 
400  ENDIF 

DO  450  J4=1,V 
HH1=B(L,J4) 

B(L,J4)=B(I3,J4) 

B(I3, J4)=HH1 
450  CONTINUE 
HH1=W1(L) 

W1(L)=W1( 13) 

Wl( I3)=HH1 
500  DO  650  J5=I4,V 

T=B( J5,I3)/B(I3,I3) 

DO  600  K=I4,V 
B(  J5  ,K)=B(  J5,K)-TlVB(I3,K) 

600  CONTINUE 

Wl( J5)=W1( J5)-T*W1(I3) 

650  CONTINUE 
655  CONTINUE 

IF  (B(V,V)  .EQ.  0.0)  THEN 
GO  TO  850 

ELSEIF  (B(V,V)  .NE.  0.0)  THEN 
GO  TO  660 
660  ENDIF 

D(V)=W1(V)/B(V,V) 

DO  680  I5=(V-1) ,1,-1 
S6=0.  0 

DO  670  J6=I5+1,V 
S6=S6+B(I5,J6)*D(J6) 

670  CONTINUE 

D(I5)=(W1(I5)-S6)/B(I5,I5) 

QlrirfrMrit'tririt 

C  WRITE(2,675)  D(I5),I5,D(V) 

C  675  FORHAT('D=' ,F25.  18,2X,I5,F25.  18) 
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680  CONTINUE 

**Calculation  of  Electric  Capacitance  using  Eq.  (57)** 
CAPE=0.  0 

DO  750  M3=1,MMAX 
DO  700  N3=1,NMAX 
N4=2*N3-1 
Q 1=MMAX*( M3 - 1 ) +N3 
KQ=(KM(M3)**2+KN(N4)**2)**0.  5 
Kl=COSH(KQ*Hl) 

K2=C0SH(KQ*H2) 

K3=C0SH( KQ*H3 ) 

S1=SINH(KQ*H1) 

S2=SINH(KQ*H2) 

S3=SINH(KQ*H3) 

S13=SINH(KQ*(H1+H3)) 

C1E=K1*S2*K3+ER*K2*S13+(ER**2)*S1*S2*S3 
C2E=D(Q1)*KQ/(KM(M3)*KN(N4)) 

C3E=K1*S2+ER*K2*S 1 

C4E=SIN(KM(M3)*A/2)*SIN(KN(N4)*Dl/2) 
NCAPE=C1E*C2E*C4E/C3E 
CAPE=CAPE+NCAPE 
C  WRITE(2,690)  CAPE,NCAPE 

C  690  F0RMAT(2( 1X,F30. 25)) 

CAPE=CAPE/ ( PI*9*10**9 ) 

700  CONTINUE 
750  CONTINUE 

WRITE(2,910)  CAPE 

910  FORMAT( 'ELECTRIC. CAPACITANCE=' ,F40, 35) 

C 

C  **Calculation  of  Magnetic  Capacitance  using  Eq. (57)** 
CAPM=0. 0 

DO  800  M5=1,MMAX 
DO  770  N6=1,NMAX 
Q2=MMAX*( M5 - 1 )+N6 
N5=2*N6 

KQ=( KM( M5 ) **2+KN ( N5 ) **2 ) **0 .  5 
Kl=COSH( KQ*H 1 ) 

K2=COSH(KQ*H2) 

K3=COSH(KQ*H3) 

S1=SINH(KQ*H1) 

S2=SINH(KQ*H2) 

S3=SINH(KQ*H3) 

S13=SINH(KQ*(H1+H3)) 

C 1M=K1*S2*K3+ER*K2*S 1 3+( ER**2 ) *S 1*S  2*S3 
C2M=D(Q2)*KQ/(KM(M5)*KN(N5)) 

C3M=K1*S2+ER*K2*S 1 

C4M=SIN(KM(M5)*A/2)*SIN(KN(N5)*Dl/2) 
NCAPM=C1M*C2M*C4M/C3M 
CAPM=CAPM+NCAPM 
CAPM=CAPM/ ( PI*9*10**9 ) 

C  WRITE(2,765)  CAPM.NCAPM 

C  765  FORMAT(2(1X,F30.  25)) 

770  CONTINUE 
800  CONTINUE 

WRITE(2,810)  CAPM 
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810  FORMATC 'MAGNETIC.  CAPACITANCE=' ,F40. 35) 
**Calculation  of  Gap  Capacitance** 

CGAP=( CAPE - C APM ) / 4 

**Calculation  of  Paracitic  Capacitance** 

CPAR=(CAPM-(Dl-G)*C0)/2 

WRITE(2.820)  CGAP 

820  F0RMAT( 'GAP1 ,1X, 'CAP=' ,F40.  35) 
WRITE(2,830)  CPAR 
830  FORMATC 'PAR. CAP=' ,F40. 35) 

850  WRITE(2,900) 

900  FORMATC IX,' MATRIX  CINGULAR  NO  SOLUTION') 
STOP 
END 
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